Anomalous scaling and spin-charge separation in coupled chains 
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We use a bosonization approach to show that the three dimensional Coulomb interaction in coupled 
metallic chains leads to a Luttinger liquid for vanishing inter-chain hopping t±, and to a Fermi liquid 
for any finite t±. However, for small t j_ ^ the Greens-function satisfies a homogeneity relation with 
a non-trivial exponent j c b in a large intermediate regime. Our results offer a simple explanation 
for the large values of j c b inferred from recent photoemission data from quasi one-dimensional 
conductors and might have some relevance for the understanding of the unusual properties of the 
high-temperature superconductors. 
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In a recent Letter [1] Chakravarty and Anderson poin- 
ted out that their interlayer-tunneling theory of high- 
temperature superconductivity [2] is based on the fun- 
damental assumption that the single-particle Greens- 
function G(k,u>) in a sufficiently large regime of wave- 
vectors k and frequencies u> satisfies the homogeneity re- 



lation G(sk, su) 
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G(k, w), with an exponent j > 0. 



It is well known that the Greens-function of a Luttinger 
liquid (LL) satisfies such a scaling law. In a Fermi li- 
quid (FL), however, the anomalous exponent j vanishes. 
There exists now general agreement that for regular in- 
teractions the LL exists only in dimensions d = 1. To 
stabilize the LL in higher dimensions, it is necessary to 
invoke singular interactions [3] , the existence of which in 
realistic experimental systems remains controversial [4]. 
A generally accepted microscopic justification for the as- 
sumption of anomalous scaling made by Chakravarty and 
Anderson does not exist. In this Letter we offer a possible 
justification. 

We proceed in three steps: We first consider a highly 
anisotropic metal consisting of an array of chains embed- 
ded in a three dimensional lattice. Assuming that the 
electrons interact via Coulomb forces, we show that for 
vanishing inter-chain hopping t±_ the system is a LL, and 
calculate the anomalous exponent j c i,. At parameters re- 
levant in recent photoemission experiments [5], we find 
that j c i, is of the order of unity. As a second step, we 
consider the same system with t±_ ^ 0, assuming that t±_ 
is small compared with the Fermi energy Ep. We show 
that, in spite of the fact that for any finite t±_ the system 
is a FL, the Greens-function for small enough \tj_/Ep\ ex- 
hibits anomalous scaling in a large intermediate regime, 
with the same exponent j c i, as for t±_ = 0. Hence, for 
sufficiently small inter-chain hopping the experiments in- 
deed measure the properties of the LL-fixed point. In this 
sense t±_ is an irrelevant perturbation. We argue that 
our result is to a large extent model independent, and 
that there exists at least a whole class of Fermi liquids 
with small quasi-particle residue Z < 1 which exhibits 



the characteristic LL-features in a large intermediate re- 
gime of momenta and frequencies. To strengthen this 
hypothesis, we finally consider an exactly solvable one- 
dimensional FL with exponentially small Z, and show 
that the spectral function of this model is practically in- 
distinguishable from that of a LL. 

We start with an anisotropic 3c? metal consisting of a 
quadratic array of Id chains where a±_ is the distance 
between the chains. As the fundamental force between 
electrons in a metal is determined by the Coulomb po- 
tential e 2 /|r — r'|, we do not replace the bare interaction 
by hand by an effective short-range interaction, but solve 
the screening problem explicitly. For a single chain the 
one-dimensional Fourier transform U(q\\) requires a re- 
gularization at short distances, which will be discussed 
below. For small one obtains a logarithmic singularity 
in U (</||) such that the anomalous dimension diverges and 
the system is not a LL [6]. The Coulomb coupling bet- 
ween the chains removes this singularity and the three 
dimensional system without hopping between the chains 
is a LL [7]. 

For the individual chains in the high density limit we 
can ignore scattering processes with large momentum 
transfer and linearize the energy dispersion around the 
two Fermi points. Then without hopping between the 
chains the electron number in each chain is conserved 
separately and we have a problem of coupled Tomonaga- 
Luttinger models (TLM) which can e.g. be solved with 
the help of a canonical transformation [8]. As we later 
address also the problem including transverse hopping 
we use here the results of the bosonization approach for 
dimensions d > 1 [9-11]. The expression for the one- 
particle Greens-function of a given chain involves the 
RPA screened interaction 
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where the parenthesis < ... > denotes the average 
over the two-dimensional first Brillouin zone (BZ), i.e. 
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<•••>= (&) 2 n;„\Au.... Here n (9||,z) = 
(2hvp q 2 /ir)/[(hvp q\\) 2 — z 2 ] is the Id non-interacting den- 
sity response-function, where vp is the Fermi velocity. 
U{q\\, qj_) is the mixed 2d-lattice and lrf-continuous Fou- 
rier transform of the Coulomb potential, 

f „. e -«'[3||»+q-L r J-] 

The rj_ = term we regularize by assuming a Gaussian 
form of the one electron wave function perpendicular to 
the chains with a decay length do [12]. As the effective 
interaction involves the Id density response function, the 
one-particle Greens-function has a form very similar to 
the single chain model [13,14]. If a = +(— ) denotes right 
(left) moving electrons, the part of the Greens-function 
describing photoemission can be written as G^(x,t) = 
[G<0M)] (O) exp[Q<0M)L where [G<(M)] (0) is the non- 
interacting Greens-function, and Q<(x,t) = 5^(0,0) — 
S< (x,t), with 

S<(x,t) = -~ V — \ e - ia n* f< e ->(q)3||* > 
2 L z — ' on I L 
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-e*'^*"*] + cos( g | F ) < 2s 2 (q)e 8 '^( q )«ll t >} , (3) 

where L is the length of the chains, vp(q) = 
v F [l + 2U(q\\,qi_)/(irhv F )]i, and s 2 (q) = [v F (q) - 
v F ] 2 /[iv F (q)v F ]. The only difference to the single-chain 
result is the average < . . . > over the BZ. The relevant 
dimensionless coupling constant follows from Eqs. (2) 
and (3) as g = e 2 /(irhv F ). The anomalous dimension 
j cb is given by j cb = ilim 3|| _ < 2s 2 (q) >, and is di- 
verging for a single chain with Coulomb forces [6]. In 
order to discuss the result for coupled chains we have 
to evaluate U{q\\ , q_i_). For arbitrary qj_ this has to be 
done numerically where we use an Ewald summation 
technique [15]. For |qj_|aj_ <C 1, we may approximate 
U(q\\ , qj_) ~ 47re 2 /(aj_q) 2 , the usual 3c? Fourier trans- 
form of the Coulomb potential. As mentioned before, the 
rj_ = contribution in Eq. (2) has a term — cln((/||ao), 
diverging for — ► 0. This divergence is exactly cancelled 
by a term cln((/||aj_) from the rj_ ^ part of the sum, 
leaving a contribution cln(aj_/ao). For every finite value 
of |qj_| therefore s 2 (q) tends to a finite value for — ► 
and the system is a LL. 

The momentum integrated spectral function p(u>), 
which apart from a one-electron dipole matrix element 
determines angular integrated photoemission, is alge- 
braically suppressed near the chemical potential p(u>) oc 
|w| 7cb . Recent photoemission studies of quasi-one dimen- 
sional conductors suggest values for the anomalous di- 
mension in the range 1.0 ±0.2 [5]. This would be hard to 
reconcile with a model involving short range interactions, 
as e.g. the anomalous dimension in the one dimensional 
Hubbard model never exceeds 1/8. Our treatment using 



real Coulomb forces shows that this quite easily leads to 
anomalous dimensions in the experimental range, when 
the value of the dimensionless coupling g is estimated 
for these systems [16]. In Fig. 1 we show results for 
j c i, as a function of g for two values of a±/ao using the 
numerical evaluation of the sum in Eq. (2). As a com- 
parison we also show the result when U{q\\, qj_) is repla- 
ced by 47re 2 /(aj_q) 2 , which provides the correct descrip- 
tion for g\ ln(aj_/ao)| <€i 1. In this limit the averages 
over the BZ can be evaluated analytically for t = [17]. 
This shows explicitly that the BZ averaging in Eq. (3) 
leads to an ultraviolet cutoff for the </||-summation which 
in this limit is given by the Thomas-Fermi wave-vector 
k = a^lSirg]? = a~^[%e 2 / (Tiv F )] ~ and can be determi- 
ned numerically in the more realistic parameter regimes. 
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FIG. 1. The anomalous dimension j c b as a function of the 
coupling constant for two values of a±/ao (dashed and dotted 
lines) and f7(g||,qj_) = 47re 2 /(aj_q) 2 (solid line). 

In order to evaluate the momentum dependent spectral 
function p < (gy , u>) we have used the numerical techniques 
developed in Ref. [14]. It is not obvious from Eq. (3) if a 
charge peak occurs as the charge- velocity v F (q) depends 
on qj_ and is subject to the BZ averaging. Because the 
spin velocity v F is independent of qj_ this leads for j c i, < 
i to a sharp threshold singularity at ui = — v F \q\\\, as in 
the lrfTLM [14,18]. For j c i, > ^ the singularity is washed 
out, but the threshold survives. To clearly demonstrate 
spin-charge separation, we therefore present in Fig. 2 
spectra for j c i, = 0.407, corresponding to g = 1. 

Experimentally the inter-chain hopping t±_ can never 
be completely turned off. In the renormalization group 
sense, t±_ is a relevant perturbation and leads to a fi- 
nite quasi-particle residue Z > [19]. We now use the 
bosonization approach in d > 1 [9-11] to show that for 
small 6 = \ti_l E F \ the behavior of the Greens-function 
at experimentally measurable length scales is completely 
dominated by the LL-exponent. Thus, from the experi- 
mental point of view, t±_ is an irrelevant perturbation. A 
crucial step in the bosonization approach is the subdi- 
vision of the Fermi surface into a set of disjoint patches 
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FIG. 2. The Lorentzian broadened spectral function 
for the coupled chains as a function of lo for different gy . 
The parameters are g = 1, a±/ao = 3 and the broadening 
x/(vfk) = 0.05. u; is measured relative to p and gy relativ to 
kF- 

For tj_ = these patches can be identified 
with the two sheets of the Fermi surface, K± = 
\k x = ±kp ; IfcjJ, \k z \ < 7r/aj_}, where we have assumed 
that the chains lie in the ^-direction. It is not necessary 
to divide K± into smaller patches, because the velocities 
v„ = Vkfklkgx = ihvp\i x are constant on K±. Here 
\i x is a unit vector in the ^-direction. Thus, Eq. (3) 
can be obtained as a special case of the general bosoniza- 
tion formula for the Greens-function in d > 1, which for 
arbitrary patch topology has been written down in Ref. 
[11]. In contast to the tj_ = calculation, this approach is 
approximate in character. In the quasi one-dimensional 
materials discussed in Ref. [16] the intra-chain hopping 
t x is an order of magnitude larger than the inter-chain 
hopping t y = which is again a factor of ten larger 
than t z . We therefore assume transverse hopping only 
in the ^-direction. In this case the sheets K± acquire 
a small modulation along the fc^-direction, with ampli- 
tude proportional to 6k p. To describe this modulation 
in terms of patches, let us assume that the band struc- 
ture is such that the distortion of patches K± can be 
approximated by wedges with hight 6k p and peak on 
the fc^-axis. The Fermi surface can then be described 
by four patches K a , a = 1, . . . , 4, on which the velocities 
v„ = [±kj; cos 6±\i y sin 6]hvp are constant. Using the ge- 
neral bosonization formula for the Greens-function in an 
arbitrary patch-geometry [11], we find after some lengthy 
but straight-forward algebra [17] that for < 6 <C 1 



R a = lim 5<(0,0) 
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and for 1 < k\r»\ < 6' 1 
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let cutoff, which for small g \ ln(aj_/ao)| can be identified 
with the Thomas- Fermi wave- vector k, and otherwise is 
in the range between aj 1 and a^ 1 . The important point 
is now that the prefactor of the logarithmic terms in Eqs. 
(4) and (5) is given by the anomalous dimension j c i, of 
the LL that would exist for 9 = 0. Because R a is fi- 
nite, the system is a FL [11], with quasi-particle residue 
Z = e Ra oc 6~ fch . Thus, Z vanishes as a non-universal 
power-law as 6 — ► 0, with the exponent given by the 
anomalous dimension j c i, of the corresponding LL. The 
transition LL — ► FL is due to the symmetry-breaking 
of the shape of the Fermi surface [9], which induces a 
coupling of the phase space associated with the q x - and 
^-integrations. Combining Eqs. (4) and (5), we see that 
< 5«( r lh 0) = — Tc6 [ln(/c | ) + ci -c 2 ]. This is precisely the 
characteristic LL-behavior, so that for 1 <C k\r\\ \ <C 
the equal time Greens-function is of the LL-form and sa- 
tisfies the scaling law G a (r/s,0) = s 3+ " tcb G a (r, 0), with 
an exponent j c i, that is determined by the corresponding 
LL at the same value of the coupling constant. 

The finite-t behavior of S< , t) cannot be obtained 
analytically. However, we expect that anomalous scaling 
holds also at finite t, and that the appearance of charac- 
teristic LL-features in our system is a model-independent 
property of a whole class of Fermi liquids with small 
quasi-particle residue. A detailed numerical study of this 
point for our anisotropic Coulomb-model with t±_ ^ 
will be presented in Ref. [17]. 

To give further evidence in favour of the generality of 
our observation, we now briefly discuss a very different 
FL with ^<1, and show that it has similar properties 
as our anisotropic Coulomb-model. Suppose we replace 
in the TLM the usual dimensionless couplings g 2 and g^ 
by F q = g{q/K) n e"«/ K , with < r) < 1, g > 0, and 
k <C kp. The result for G^(x,t) is of the same form as in 
Eq. (3), except that the BZ averages should be omitted, 
and vp(q) should be replaced by vp(q) = vp[l + 2F q ]i . 
While for rj = the q — ► limit of F q is finite and leads 
to a logarithmic divergence of 5^(0,0) as L — ► oo and 
LL behavior, this divergence is suppressed for rj > and 
one obtains a FL with quasi-particle residue Z = e Ra . 
Obviously, this FL is obtained from the TLM by defor- 
ming the potential energy, whereas the distortion of the 
shape of the Fermi surface in our three dimensional mo- 
del discussed above is a kinetic energy effect. To examine 
if the Greens-function satisfies anomalous scaling in some 
parameter regime, we calculate 

G<(x/ S ,t/ S ) = S [G<(x,t)]^Ze- s ^^ s ^^ . (6) 

We find for the quasi-particle residue Z = e~^ g ^l^^ g \ 
and for S^(x, 0) in the large distance limit \kx\ 1 



S<(x,0) 
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where c\ and c 2 are numerical constants of the order of 
unity, and r*y 



rl + * 2 ~ 2 



r-y [20]. Here k is an ultravio- 



f]{g)\K,x\^Ji) ' 
where to leading order in rj the exponent j(g) can be 
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identified with the anomalous dimension of the corre- 
sponding TLM with coupling constants g 2 = g^ = g, 
i.e. j(g) = i[(l + 2g)i - 1] 2 /(1 + 2g)i . The exponent 
rj(g) is a featureless monotonically decreasing function of 
g, varying between fj(0) = 2rj and 77(00) = rj/2. Note 
that for any g > and sufficiently small rj the quasi- 
particle residue becomes exponentially small. It is now 
trivial to prove the existence of anomalous scaling in an 
exponentially large regime. Writing | Ka; | — ^ = e ~ f i iD \ KX \ 
and assuming | In | /car 1 1 <C 1, Eq. (7) implies to lea- 
ding order S< (x/s, 0) k> S< (x, 0) — j(g) In s. Using then 
Eq. (6), we conclude that G<(x/s,0) « s 1+ f^ ~>G< (x, 0) 
in the exponentially large regime 1 <C \kx\ <C e 1 /^ 3 ). 
The restriction due to the lower limit holds also in a 
LL, but in the FL anomalous scaling breaks down at 
x max (g) = K~ 1 e 1 l ft ( g \ Obviously the anomalous scaling 
regime exists only for rj <C 1. Note also that at strong 
coupling x max (g) is much larger than at weak coupling, 
in agreement with the intuitive picture that strong corre- 
lations should lead to more pronounced anomalous pro- 
perties. 
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FIG. 3. The Lorentzian broadened spectral function for 
the Id FL as a function of lo for different q. The parameters 
are: r] = 0.1, g = 2 and x/("f K ) = 0.05. The small 
oscillations are a finite-size effect. 

At finite t we have used again the technique of Ref. 
[14] to calculate the spectral function. In Fig. 3, we 
present the spectral function of the FL for rj = 0.1, which 
is to the accuracy of the drawing indistinguishable from 
the spectral function of a TLM with g 2 = #4 = g [21]. 
Because we know that the spectral function of the TLM 
satisfies anomalous scaling, we conclude for rj > the 
FL must have this property as well. From Fig. 3 it is 
obvious that our FL exhibits spin-charge separation in 
the scaling regime. 

Our result provides a very simple microscopic justi- 
fication for the assumption of anomalous scaling made 
by Chakravarty and Anderson [1] to explain the unusual 
properties of the high-temperature superconductors. As- 
suming that due to strong correlations and anisotropy 



these systems are Fermi liquids with small quasi-particle 
residue, we predict that they should exhibit all the cha- 
racteristic LL-features in the experimentally accessible 
range of wave- vectors and frequencies. Thus, although 
dimensionality and transverse hopping are relevant in the 
renormalization group sense, they can be considered as 
irrelevant as far as the experiments are concerned. 

One of us (K.S.) would like to thank D. Baeriswyl and 
G. M. Eliashberg for stimulating discussions. 



[1] S. Chakravarty and P. W. Anderson, Phys. Rev. Lett. 

72, 3859 (1994). 
[2] S. Chakravarty et al, Science 261, 337 (1993). 
[3] P. W. Anderson, Phys. Rev. Lett. 64, 1839 (1990); 65, 

2306 (1990); 66, 3226 (1991). 
[4] J. R. Engelbrecht and M. Randeria, Phys. Rev. Lett. 65, 

1032 (1990); 66, 3225 (1991). 
[5] B. Dardel et al, Europhys. Lett. 24, 687 (1993); M. Na- 

kamura et al, Phys. Rev. B 49, 16191 (1994); R. Claes- 

sen, private communication. 
[6] H. J. Schulz, Phys. Rev. Lett. 71, 1864 (1993). 
[7] H. J. Schulz, J. Phys. C 16, 6769 (1983). 
[8] K. Penc and J. Solyom, Phys. Rev. B 47, 6273 (1993). 
[9] F. D. M. Haldane, Helv. Phys. Acta. 65, 152 (1992); in 

Proceedings of the Int. School of Physics "Enrico Fermi" , 

Course 121, (North Holland, Amsterdam, 1994). 
[10] A. Houghton and J. B. Marston, Phys. Rev. B 48, 7790 

(1993) ; A. Houghton et al, ibid. 50, 1351 (1994); A. H. 
Castro Neto and E. Fradkin, Phys. Rev. Lett. 72, 1393 

(1994) ; Phys. Rev. B 49, 10877 (1994); and preprint. 
[11] P. Kopietz and K. Schonhammer, preprint, July 1994. 
[12] This amounts to -> (■Ka 2 )~ 1 f (ixj_e" x i /a o(x 2 + 

xi)" 1/2 . 

[13] F. D. M. Haldane, J. Phys. C 14, 2585 (1981). 

[14] V. Meden and K. Schonhammer, Phys. Rev. B 46, 15753 

(1992) ; K. Schonhammer and V. Meden, ibid. 47, 16205 

(1993) . 

[15] J. M. Ziman, Principles of the Theory of Solids, (Cam- 
bridge University Press, Cambridge, 1964), p. 39. 

[16] D. Jerome, Science 252, 1509 (1991). 

[17] P. Kopietz, V. Meden and K. Schonhammer, to be pu- 
blished. 

[18] J. Voit, Phys. Rev. B 47, 6740 (1993). 

[19] C. Castellani et al, Phys. Rev. Lett. 69, 1703 (1992). 

[20] [G<(r,i)] (0) is proportional to <5 (d_1) (r a j_), where r ai is 

the d — 1 dimensional projection of r perpendicular to 

v„ . The r-dependence of Q a enters therefore only via 

the component of r parallel to v a . 
[21] For 7/ = 0.1 we found numerically that 7 is given more 

accurately if one takes the maximum value of F q , i.e. 

n = gi = gv v e~ v - 



4 



